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I.  INTRODUCTION 


It  Is  well  known  that  radiation  from  a  nuclear  explosion  can  generate  a 

thin  layer  of  heated  air  adjacent  to  the  ground.^  The  interaction  of  the 

nuclear  blast  wave  with  this  thermal  layer  results  in  a  precursor  wave  system, 

which  affects  dust  cloud  generation  and  the  aerodynamic  load  on  ground 

structures.^  A  general  discussion  of  this  interaction  has  been  provided  by 

Hess,^  Experimental  studies  of  the  interaction  of  weak  (M  <  1.14)  moving 

s 

shocks  with  thin  thermal  layers  have  been  reported  in  Refs.  3  and  4. 

Detailed  numerical  code  calculations  of  the  inviscid  interaction  between 
an  incident  shock  and  a  thin  thermal  layer  of  semi -infinite  streamwise  extent 
have  recently  been  reported.^  The  present  study  is  an  attempt  to  deduce 
properties  of  the  latter  interaction  by  simple  analytic  methods  and  to  compare 
with  the  numerical  results  of  Ref.  5. 
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II.  THEORY 


Consider  a  normal  shock  wave  moving  over  a  thin  thermal  layer  of  semi- 
infinite  extent  and  assume  ideal  gases.  The  Initial  flow  field,  in  laboratory 
stationary  coordinates,  is  Illustrated  in  Fig.  1.  For  generality,  the  ratio 
of  specific  heats  y  is  allowed  to  differ  in  the  thermal  layer  and  in  the 
external  flow.  After  the  initial  transient,  the  interaction  can  be  charac¬ 
terized  as  either  unseparated  or  separated.  These  flows  are  steady  and 
unsteady,  respectively,  in  Incident  shock  stationary  coordinates.  Both  cases 
are  considered  herein.  Flow  velocity  and  Mach  number  are  denoted  by  u  and  M, 
respectively,  in  the  laboratory  coordinate  system  (Fig.  1),  and  by  u  and  M  in 
the  incident  shock  stationary  coordinate  system  (Fig.  2).  Note  that 

”g  =  “s^®l  “ 

A.  UNSEPARATED  FLOW  REGIME 


We  consider  unseparated  steady  flow  in  a  shock  stationary  coordinate 
system  (Fig.  2).  The  static  pressure  downstream  of  the  incident  shock  is^ 


P2  ZYjMJ  -  (Yj  -  1) 

^  Tpn. 


(1) 


Is^ 


The  stagnation  pressure  in  the  downstream  portion  of  the  thermal  layer 
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As  noted  in  Ref.  2,  the  interaction  between  the  shock  and  thermal  layer  will 
remain  unseparated  and  steady  (in  shock  fixed  coordinates)  provided 
P5  j./p4  >  P2/Pi«  The  flow  becomes  separated  and  unsteady  when  < 


'  ^ 


P2/p^*  The  boundary  between  the  separated  and  unseparated  flow  regimes  Is 
obtained  by  equating  Eqs.  1  and  2  and  noting  =  (aj^/a^)Mp  The  results  are 
given  In  Fig.  3  and  agree  with  similar  results  In  Ref.  2.  The  unseparated 
flow  regime  occurs  for  values  of  a^/a^  and  Mj  near  one.  Increased  values  of 
a^/a^  and/or  lead  to  a  separated  flow,  which  is  discussed  in  the  next 
section. 

B.  SEPARATED  FLOW  REGIME 

Numerical  results  from  Ref.  3  and  an  analytic  model  are  presented  herein. 
1 .  NUMERICAL  RESULTS 

Detailed  numerical  calculations  of  the  interaction  between  an  incident 
shock  and  a  thermal  layer  have  been  presented  In  Ref.  5  for  various  Incident 
shock  Mach  numbers  and  speed  of  sound  ratios  a^/a^.  Specific  heat  ratio 
values  »  1.4  and  *  5/3  (approximately)  were  assumed  so  that  the 
interface  between  the  fluid  In  regions  1  and  4  was  well  defined.  Viscous 
effects  were  neglected. 

Results  from  Ref.  5  are  given  in  Table  1  and  Fig.  4.  Figs.  4(a)  and  4(b) 
provide  density  contours  at  a  fixed  Instant  of  time  (Xg  =*  20.1  cm),  and  Fig. 
4(c)  provides  the  corresponding  pressure  distribution  at  three  heights  above 
the  ground.  Initial  conditions  were  M^  *  2.61  and  a^/aj^  «  2.26.  The  latter 
w?'  termed  Case  II  In  Ref.  5.  A  schematic  representation  of  the  flow  field, 
inferred  from  Fig.  4,  Is  given  In  Fig.  5(a)  using  coordinates  wherein  the 
incident  shock  is  stationary.  Shock-heated  gas  originally  from  region  1  moves 
forward  of  the  Incident  shock  and  separates  the  thermal  layer  gas  flow  from 
the  wall.  The  Interface  between  gas  from  region  1  and  gas  from  region  4  Is 
Indicated  by  a  dashed  line  In  Fig.  5(a).  The  Interface  location  at  the  wall 
Is  denoted  x^.  The  quantities  x^p,  Xg,  and  Xg^,  in  Fig.  5(a),  denote  the 
streamwise  location  of  a  stagnation  point,  the  incident  shock,  and  the  leading 
edge  of  the  normal  shock  in  the  thermal  layer,  respectively.  The.  correspond¬ 
ing  locations  of  these  stations  are  included  in  Fig.  4.  The  variation  of 
these  locations,  with  time,  deduced  from  the  data  reported  in  Ref.  5,  are 
given  in  Fig.  6.  The  local  slope  of  the  curves  in  Fig.  6  define  local 
velocity.  These  results  indicate  that,  after  an  initial  transient,  the 
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Fig,  4a. 


Numerical  Results  of  Ref.  5  for  Shock- Thermal  Layer  Interaction 
Case  II,  Xg  -  20.1  cm.  Initial  conditions  noted  in  Table  1. 
(Contours  of  constant  density,  overall). 
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Fig.  4b.  Numerical  Results  of  Ref.  5  for  Shock-Thermal  Layer  Interaction 
Case  II,  Xg  ■  20.1  cm.  Initial  conditions  noted  in  Table  1. 
(Contours  of  constant  density,  detail). 
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Flg«  4c.  Numerical  Results  of  Ref.  5  for  Shock- Thermal  Layer  Interaction 
Case  II,  Xg  ■  20.1  cm.  Initial  conditions  noted  in  Table  1. 
(Overpressure  at  three  heights  (pj^  ■  1.013  x  10®  dynes/cm^)). 


Interface  Trajectories  Based  on  Numerical  Solution  of  Kef 
Initial  conditions  noted  In  Table  1.  (Case  II  (aA/ai  =  2 


Interface  Trajectories  Based  on  Numerical  Solution  of  Ref.  3 
Initial  conditions  noted  In  Table  1.  (Case  VII  (aASt  =  2.92 


Interface  location  moves  with  constant  speed.  The  wall  stagnation  point 
location  moves  with  a  velocity  approximately  equal  to  that  of  the  Incident 
shock  Xg  and,  hence.  Is  stationary  In  Incident  shock  coordinates.  Similarly, 
the  leading  edge  of  the  normal  shock  in  the  thermal  layer  Xg^  moves  with 
approximately  the  same  velocity  as  x^.  Thus,  the  separation  distance  Xg^  -  Xj^ 
remains  constant.  The  flow  In  a  coordinate  system  wherein  Xg^  and  x^^  are 
stationary  Is  Illustrated  In  Fig.  5(b).  The  flow  In  region  5  Is  steady  In 
this  coordinate  system.  These  results  suggest  an  analytic  model  described  In 
the  next  section. 

2.  ANALYTIC  MODEL 

We  consider  the  flow  In  an  incident  shock  stationary  coordinate  system  as 
shown  In  Fig.  5(a).  The  stagnation  point  Xgp  Is  assumed  to  be  stationary,  and 
the  separation  distance  Xg^-x^^  is  assumed  to  remain  constant  for  cases  when  a 
shock  develops  In  the  thermal  layer.  Flow  conditions  In  each  of  the  regions 
noted  In  Fig.  5  are  obtained  from  the  following  expressions. 

a.  Region  2 

Region  2  Is  the  uniform  region  downstream  of  the  incident  shock.  Normal 
shock  relations  Indicate^ 


P2  2TjM2  -  (Y^  -  1) 
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where  P2  ^  Is  the  stagnation  pressure  In  region  2.  Since  Xgp  Is  stationary  In 
the  present  coordinate  system. 


Psp'^Pl  “  P2,t/Pl 


(3c) 


The  possibility  that  the  flow  at  p^p  was  processed  by  the  oblique-nonnal  shock 


structure  at  the  base  of  the  Incident  shock  is  ignored  in  Eq,  3(c).  However, 


the  effect  on  p^p  should  not  be  large. 


b •  Region  3 


Conditions  in  region  3  are  assumed  to  result  from  a  steady  isentropic 
expansion  from  conditions  in  region  2.  Expressions  of  Interest  are 
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where  Q  denotes  dynamic  pressure  in  laboratory  coordinates. 


Region  5 


For  cases  in  which  a  shock  develops  in  the  thermal  layer,  the  separation 
distance  x^^-x^  is  assumed  to  remain  constant  (Fig.  6).  The  flow  in  region  5 
is  then  equivalent  to  steady  supersonic  flow  with  Mach  number 


- 


M^  +  (uj./aj) 


St  a. 


over  a  blunt  body  defined  by  the  interface  at  x^.  For  cases  in  which  a  shock 
does  not  develop,  the  flow  in  region  5  is  equivalent  to  subsonic  flow  at  Mach 
number  over  a  blunt  body  defined  by  the  -interface.  The  pressure  pj^ 

corresponds  to  the  stagnation  pressure  at  the  nose  of  the  equivalent  body  and 
is  found  from 
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where  and  u^^/a^  can  be  used  interchangeably  (Eq.  5). 
d.  Method  of  Solution 

The  interaction  is  uniquely  determined  by  specification  of  the  initial 
conditions  M^f  Y^,  a^/ap  and  h.  (The  quantity  h  is  not  needed,  however, 
for  the  present  model.)  Estimates  for  flow  properties  can  be  obtained  from 
the  previous  equations  if  the  boundary  conditions  at  the  interface  location  x. 
are  specified.  In  the  present  model,  we  assume  that  conditions  in  regions  3 
and  5,  at  Xj^,  are  related  by 

Pi  “  P3  (7a 


The  arbitrary  constant  k  is  Introduced  to  permit  matching  of  the  analytic 
model  with  numerical  code  calculations.  A  value  of  k  <  1  can  be  Interpreted 
as  indicating  a  transverse  velocity  in  region  3  at  the  interface  (Fig.  7).  A 
simple  procedure  for  evaluating  the  previous  equations  is  to  specify  y  p  Y^» 
Ml  =  M^  ,  and  M3.  M^^.  (or,  equivalent,  u^/a^)  is  found,  by  iteration,  by 


Boundary  Conditions  Near  Interface  In  Incident  Shock  Stationary 
Coordinate  System,  k  =  cos  9. 


equating  Eqs.  6  and  4(a).  The  corresponding  values  of  and  are 

obtained  from  Eqs.  5  and  7(b),  namely 


aj/*! 


[(M  +  (ku-/a  ) ]/M 

1  i  1  St 


Uf/ai 


(u^/a^)(a^/a^) 


Using  the  present  procedure,  a^/a^  and  u^^/a^  are  the  only  quantities  that 
depend  on  the  value  of  k.  Numerical  results  are  given  In  Table  2  for  the  case 
k  *  1.  For  a  given  Incident  shock  M  ,  the  dynamic  pressure  ratio  Q,/Q~  has  a 
maximum  at  a  particular  value  of  a^/a^.  The  maximum  Is  characteristic  of  the 
Isentroplc  expansion  to  Values  of  In  the  range  <  1  Indicate  that 

a  shock  Is  not  formed  In  the  thermal  layer. 
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III.  DISCUSSION 


The  present  analytic  model  Is  compared  in  Table  1  vlth  numerical  data 
from  Ref.  5.  Values  of  u^/a^,  Pgp/Pl*  Pl/Pl»  from  Ref.  5,  are  given 

as  a  function  of  Incident  shock  location  Xg  for  two  Initial  conditions  (Cases 
II  and  VII).  Analytic  estimates  for  these  quantities,  based  on  k  >  1.00  and  k 
»  0.27,  are  Included  In  Table  1.  The  values  of  Pgp/Pi  In  Ref.  5  Increase  with 
Xg  and  appear  to  approach  Che  analytic  estimate.  Hence,  the  assumption 
Psp/Pl  “  P2  t^Pl  appears  reasonable.  Other  analytic  estimates  In 

Table  1  depend  on  Che  choice  for  k.  The  Initial  choice,  k  «  1.00,  led  to  an 
overestimate  of  u^/a^,  and  p^/pj.  However,  the  variation  of  these 

quantities  with  had  the  proper  trend.  The  second  value,  k  *  0.27,  was 

chosen  so  that  Che  analytic  estimate  for  pj^/p^  would  agree  with  Che  numerical 
result  In  Ref.  5  for  Case  II.  This  led  to  Improved  correlations  between  Che 
analytic  model  and  Che  Cases  II  and  VII  data  from  Ref.  5.  However,  Che 
usefulness  of  the  approximation  k  »  0.27  for  other  Initial  conditions  Is 
uncertain  and  requires  further  study. 

The  Interface  pressure  P^/p]^  was  also  computed  from  Eq.  6(a)  using  the 
value  of  M  deduced  from  Che  data  in  Ref.  5.  The  resulting  pressure  Is 

9  t 

denoted  (p^/pj)  q  In  Table  1  and  agrees  with  Che  value  of  P^/Pi  reported  in 
St 

Ref.  5.  This  calculation  confirms  Che  analytic  model  assumption  that  the 
pressure  P^/Pi  equals  the  stagnation  pressure  associated  with  (Eq.  6). 


IV.  CONCLUSION 


Numerical  solutions^  for  the  invlscld  Interaction  between  a  moving  shock 
and  a  thermal  layer  have  been  reviewed.  The  streamwlse  locations  Xgp,  x^, 
and  Xg^  were  of  particular  Interest.  It  was  determined  that  after  an  Initial 
transient  dXgp/dt  *  dXg/dt  =  This  led  to  the  assumption.  In  the  analytic 

model,  that  the  pressure  at  Xgp  can  be  equated  to  the  stagnation  pressure 
behind  the  Incident  shock  In  Incident  shock  stationary  coordinates  (Pgp  = 

P2  ^).  Similarly,  It  was  concluded  that  after  an  Initial  transient 
dXst/dt  =  dx^/dt  =  u^.  Hence,  the  pressure  pj^  at  Xj^  was  assumed.  In  the 
analytic  model,  to  be  equivalent  to  the  stagnation  pressure  on  a  blunt  body  In 
a  Mach  number  flow.  The  analytic  model  was  only  partly  successful 

In  predicting  x^  and  pj^.  The  Introduction  of  a  free  constant  k  *  u^/uj  was 
required.  This  constant  Is  related  to  the  angularity  of  the  flow  In  region 
3.  The  choice  k  «  0.27  led  to  good  agreement  between  the  analytic  model  and 
data  from  cases  II  and  VII.  However,  the  large  departure  of  k,  from  one, 
makes  the  physical  Interpretation  questionable.  Further  comparisons  with 
numerical  solutions  are  needed  to  validate  the  analytic  model. 

In  the  present  Invlscld  model,  the  scale  of  the  disturbed  flow  grows 
linearly  and  without  limit.  It  is  expected  that  viscous  effects  will  cause 
the  flow  to  approach  an  asymptotic  limit.  The  scale  and  nature  of  the  flow, 
in  this  limit,  also  requires  further  study. 
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SYMBOLS 


a  speed  of  sound 

h  initial  height  of  thermal  layer 

k  arbitrary  constant,  k  <  1 

M  Mach  number  in  incident  shock  stationary  coordinate  system,  u/a 

M  Mach  number  in  laboratory  coordinates,  u/a 

p  pressure 

Q  dynamic  pressure  in  laboratory  coordinates,  p(u)2/2 

t  time 

u,u  velocity 

x,z  streamwlse  and  vertical  directions,  respectively 

Y  ratio  of  specific  heats 

SUBSCRIPTS 

1,2, 3, 4, 5  flow  regions 

1  interface  value 

s  incident  shock 

St  shock  in  thermal  layer 

sp  stagnation  point  at  foot  of  incident  shock 

t  stagnation  value 

SUPERSCRIPTS 

(~)  superscript  bar  denotes  laboratory  coordinates 

(  )  absence  of  superscript  bar  denotes  incident  shock  stationary 

coordinates 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporatloa  functions  as  an  “architect-engineer"  for 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  Investigations  Is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  Individual  laboratories: 


Aerophyslcs  Laboratory;  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  Including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory;  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of-fleld-of-vlew  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  Infrared  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 


Computer  Science  Laboratory;  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spacebome 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory;  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications;  micro- 
wave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements,  diagnos¬ 
tics  and  radlometry,  mlcrowave/mllllmeter  wave  thermonlc  devices;  atomic  time 
and  frequency  standards;  antennas,  RE  systems,  electromagnetic  propagation 
phenomena,  space  communication  systems. 


Materials  Sciences  Laboratory;  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  aiuilysls  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory;  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics.  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere.  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
lost  rumenta  t Ion . 
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